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ABSTRACT 

The present work is designed to explore the evolutionary and pulsational properties of 
low-mass white dwarfs with carbon/oxygen cores. In particular, we follow the evolution 
of a 0.33-A/q white dwarf remnant in a self-consistent way with the predictions of 
nuclear burning, element diffusion and the history of the white dwarf progenitor. 
Attention is focused on the occurrence of hydrogen shell flashes induced by diffusion 
processes during cooling phases. The evolutionary stages prior to the white dwarf 
formation are also fully accounted for by computing the conservative binary evolution 
of an initially 2.5-MqPop. I star with a 1.25 MQeompanion, and period Pj= 3 days. 
Evolution is followed down to the domain of the ZZ Ceti stars on the white dwarf 
cooling branch. 

We find that chemical diffusion induces the occurrence of an additional hydrogen 
thermonuclear flash which leads to stellar models with thin hydrogen envelopes. As a 
result, a fast cooling is encountered at advanced stages of evolution. In addition, we 
explore the adiabatic pulsational properties of the resulting white dwarf models. As 
compared with their helium-core counterparts, low-mass oxygen-core white dwarfs are 
characterized by a pulsational spectrum much more featured, an aspect which could 
eventually be used for distinguishing both types of stars if low-mass white dwarfs 
were in fact found to pulsate as ZZ Ceti - type variables. Finally, we perform a non- 
adiabatic pulsational analysis on the resulting carbon/oxygen low-mass white dwarf 
models. 

Key words: stars: evolution - stars: interiors - stars: white dwarfs - stars: binaries - 
stars: oscillations 
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1 INTRODUCTION 

Low-mass white dwarf (WD) stars, having masses less than 
« 0.5 Mq, are the result of evolution of certain close bi- 
nary systems (Paczyhski 1976; Iben & Webbink 1989; Iben 
& Livio 1993). Indeed, mass-transfer episodes in binary sys- 
tems are required to form low-mass WDs. An isolated star of 
comparable low mass would take more than a Hubble time 
to evolve into a corresponding WD configuration. The low- 
mass WDs are thought to have a core composed by helium; 
their binary nature was first put on a firm observational 
basis by Marsh (1995) and Marsh, Dhillon & Duck (1995). 
Since then, low-mass WDs have been detected in numerous 
binary configurations containing usually either another WD 
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or a neutron star (see, e.g., Lundgren et al. 1996; Moran, 
Marsh & Bragaglia 1997; Orosz et al. 1999; Maxted et al. 
2000, van Kerkwijk et al. 2000). In addition, several low- 
mass WDs have been found in open and globular clusters 
(Landsman et al. 1997; Edmonds et al. 1999, 2001; Taylor 
et al. 2001). Evolutionary models for helium-core WDs have 
been calculated by Hansen & Phinney (1998), Driebe et al. 
(1998); Sarna, Ergma & Antipova (2000), Althaus, Serenelli 
& Benvenuto (2001a,b) and Serenelli et al. (2002). 

However, theoretical evidence casts some doubt on that 
all low-mass WDs be indeed helium-core WDs. In fact, Iben 
& Tutukov (1985) proposed several scenarios where low- 
mass white dwarfs could harbor cores with elements heavier 
than helium. More recently, Han, Tout & Eggleton (2000) 
have carried out new close binary evolutionary calculations 
and found that some of the presumed helium-core WDs 
in double degenerate systems may actually be WDs with 
carbon-oxygen cores. Specifically, their calculations show 
that, if the onset of mass transfer episodes takes place af- 
ter the end of central hydrogen burning during the middle 
Hertzsprung-gap stage, a carbon-oxygen WD with mass as 
low as 0.33 M Q may be formed from stable Roche lobe over- 
flow (RLOF) if the initial mass of the primary star is close 
to 2.5 Mq. 

The present work aims at exploring the evolution of 
low-mass WDs with carbon/oxygen cores in a self-consistent 
way with the predictions of nuclear burning, time-dependent 
element diffusion and the history of the white dwarf progen- 
itor. We concentrate on a carbon/oxygen-core WD remnant 
of mass 0.33 Mq . Attention is focused on the occurrence of 
hydrogen shell flashes and the role played by diffusion pro- 
cesses during cooling phases. The evolutionary stages prior 
to the white dwarf formation are also fully accounted for 
by computing the conservative binary evolution of an ini- 
tially 2.5-M Q Pop. I star with a 1.25-M0companion, and pe- 
riod Pi= 3 days. Finally, we explore the adiabatic and non- 
adiabatic pulsational properties of the resulting WD models 
and we compare them with those of helium-core WDs. 

In this work, we hope to address essentially the two fol- 
lowing questions: (i) What is the role of element diffusion 
in the evolution of low-mass, carbon/oxygen WDs?; and (ii) 
What differences in the evolutionary properties as well as 
in the pulsational signal (if these objects do actually pul- 
sate) can we expect between low-mass WDs with helium 
cores and those with carbon/oxygen cores? In Section 2, wc 
shortly describe the main physical inputs of the models. In 
section 3, we present the evolutionary results, particularly 
regarding the WD stages. Section 4 contains the discussion 
of the pulsational properties of our models. Finally, Section 5 
is devoted to summarizing our results and making some con- 
cluding remarks. 



2 INPUT PHYSICS AND COMPUTATIONAL 
DETAILS 

The evolutionary code used in this work to simulate the 
binary evolution is that described by Han et al. (2000) 
and references therein. The code uses a self-adaptive non- 
lagrangian mesh and it is based on an up-to-date physi- 
cal description such as OPAL radiative (Rogers & Iglesias 
1992) and molecular (Alexander & Ferguson 1994) opacities, 



a detailed equation of state that includes coulomb interac- 
tions and pressure ionization and nuclear reaction rates from 
Caughlan et al. (1985) and Caughlan & Fowler (1988). Con- 
vective overshooting is not considered. 

For the WD regime, we employed the LPCODE evo- 
lutionary code which has been used in previous works on 
WD evolution. The code has been developed at La Plata 
Observatory and it is described in Althaus (1996), Al- 
thaus et al. (2001a), Serenelli (2002), Althaus et al. (2003) 
and references therein. As for the constitutive physics, LP- 
CODE employs OPAL radiative opacities (including carbon 
and oxygen-rich compositions) for arbitrary metallicity from 
Iglesias & Rogers (1996). Opacities for various metallicities 
are required during the WD cooling regime in view of the 
metallicity gradients that develop as a result of gravitational 
settling (see below) . The equation of state is an updated ver- 
sion of that of Magni & Mazzitelli (1979). High-density con- 
ductive opacities and the various mechanisms of neutrinos 
emission are taken from the works of Itoh and collabora- 
tors. Hydrogen burning is taken into account by considering 
a complete network of thermonuclear reaction rates corre- 
sponding to the proton-proton chain and the CNO bi-cycle. 
Nuclear reaction rates are taken from Caughlan & Fowler 
(1988). Electron screening is treated as in Wallace, Woosley 
& Weaver (1982). 

An important aspect of the present study is the evolu- 
tion of the chemical abundance distribution caused by ele- 
ment diffusion during the whole WD stage. In our treatment 
of time-dependent diffusion we have considered the following 
nuclear species: 1 H, 3 He, 4 He, 12 C, 14 N and ie O. The chem- 
ical evolution resulting from element diffusion is described 
in LPCODE, for a given isotope i having a number density 
rii, by the continuity equation as 



dm 19/ 2 \ 



(i) 



where Wi is the diffusion velocity as given by the solution of 
the multicomponent flow equations describing gravitational 
settling, chemical and thermal diffusion (Burgers 1969). Be- 
cause we are interested in the chemical evolution occurring 
quite deep in the star, radiative levitation (which is expected 
to modify the surface composition of hot WDs) has been ne- 
glected. In terms of the gradient of ion densities, diffusion 
velocities can be written in the form 



Wi 



ions(j) 



d In n j 
dr 



(2) 



where wf stands for the velocity component due to gravi- 
tational settling and thermal diffusion and Oij represent the 
components due to the gradients in number density (chemi- 
cal diffusion). Details are given in Althaus et al. (2001a) and 
Gautschy & Althaus (2002). We want to mention that after 
computing the change of abundances by effect of diffusion, 
they are evolved according to nuclear reactions and convec- 
tive mixing. It is worth mentioning that radiative opacities 
are calculated for metallicities consistent with diffusion pre- 
dictions. In particular, the metallicity is taken as two times 
the abundances of CNO elements. 

In this work, we have followed the complete evolution 
of an initially 2.5-Mostar in a close binary system from the 
main sequence through the mass transfer episodes to the do- 
main of the ZZ Ceti stars on the WD cooling branch. Specif- 
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ically, the initial mass of the primary and secondary stars is, 
respectively, 2.5 and 1.25 Mowith an initial orbital period 
of Pi= 3 days. A solar-like initial composition (Y,Z)= (0.28, 
0.02) was adopted. The treatment of convection is that of 
the mixing-length theory with the ratio of the mixing length 
to the local pressure scale-height set to two (see Han et al. 
2000 for details). The onset of RLOF, which is assumed to 
be conservative, occurs after the central hydrogen exhaus- 
tion in the middle of the Hertzsprung-gap stage. After mass 
transfer episodes, the mass of the WD remnant is 0.33 Mq. 
As stated in the introduction, one of the aim of this investi- 
gation is the exploration of the influence that diffusion has 
on the evolution of low-mass, carbon/oxygen WDs. Thus, in 
addition to evolutionary models in which diffusion is consid- 
ered, we have computed the WD evolution for the case when 
diffusion is neglected. This allowed us to clearly identify the 
effect induced by diffusion on the evolution of such objects. 



3 EVOLUTIONARY RESULTS 

In the following, we describe the evolution on the 
Hertzsprung-Russell diagram (HRD) of the WD progenitor. 
Specifically, in FigQwe show the HRD of the primary star, 
which initially has a stellar mass of 2.5 Mq. The primary 
evolves from zero-age main sequence (point A) and burns 
hydrogen at its centre. After central hydrogen is exhausted, 
the star evolves to the Hertzsprung gap and fills its Roche 
lobe at point B, and the RLOF phase begins. The envelope 
mass is transferred to the secondary, which initially has a 
stellar mass of 1.25 Mq, till RLOF finally stops at point C. 
After the RLOF, the primary becomes a helium star (with 
a thin hydrogen layer). As a result of mass loss episodes, 
the stellar mass has decreased to M\{ — 0.3293 Mq, while 
the mass of the secondary is now Mil = 3.4207 M0and the 
orbital period has increased to Pf = 64.08 d. The helium 
star contracts and moves quickly across the HRD. Helium 
is ignited at the centre and the star experiences "breathing 
pulses" (the loops at lower-left corner of the figure; see e.g. 
Castellani et al. 1985, Han et al. 2002). After central helium- 
burning, the star has an almost pure oxygen core of 0.11 Mq. 
The evolution calculation of the resulting oxygen-core WD 
is continued from point D. 

The resulting WD binary system continues to evolve 
and the secondary will fill its Roche lobe as an asymptotic 
giant branch (AGB) star and leads to the onset of another 
RLOF. The star has a CO core and a deep convective enve- 
lope. The mass transfer is dynamically unstable and a com- 
mon envelope (CE: Paczyhski, 1976) is formed. The CE en- 
gulfs the CO core and the primary (a WD). A large amount 
of orbital energy is deposited into the CE due to the friction 
between the CE and the embedded binary. The CE may 
be ejected if the orbital energy released can overcome the 
binding energy of the envelope, and leads to the formation 
of a binary system with two WD components, i.e. a dou- 
ble degenerate. If the common envelope ejection efficiency is 
qce = 1 and thermal contribution to CE ejection is a t h = 1 
(see Han 1998; Han et al. 2002 for the definitions), the double 
degenerate system will have Mi = 0.3293 Mq, M 2 = 0.6572 
Mq, P = 0.1837 d. However, if a C E = «th = 0.75, the 
orbital period of the double degenerate will be 0.0735 d in- 
stead. 
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Figure 1. Hertzsprung-Russell diagram for the primary star dur- 
ing the binary evolution stage and core helium burning stages. 
Labels along the track denote specific stages during the evolution 
and are discussed in the text. 

In the following, we describe the main evolutionary re- 
sults for the WD evolution of the 0.33-Mooxygen-core rem- 
nant. In Figs.|2]and|j|]we show the HRD for the WD remnant 
when element diffusion is neglected and when diffusion is 
considered, respectively. Labels along the track from A to H 
represent evolutionary stages for models with diffusion, the 
main characteristics of which are listed in Table 1. In both 
cases, evolution is characterized by the occurrence of sev- 
eral episodes of thermonuclear flashes during which the star 
experiences extensive loops in the Hertzsprung-Russell dia- 
gram. Each of these loops is due to unstable nuclear burning 
at the bottom of the hydrogen-rich envelope. Note that when 
diffusion is considered, the model experiences an additional 
thermonuclear flash, as compared with the situation of no 
diffusion. This last flash is triggered by chemical diffusion 
which, as the WD evolves along the cooling branch, car- 
ries some hydrogen inwards to hotter layers, where it burns 
unstably 1 . As we shall see, the occurrence of this diffusion- 
induced thermonuclear flash is critical regarding the subse- 
quent evolution of the remnant even during its final cooling 
phase. It is worth noting that the diffusion-induced flash is 
the most energetic one and as a result the star, for a brief 
time interval, reaches giant dimensions again (point D). Af- 
ter this flash episode, the remnant eventually settles on its 
terminal cooling track. 

The time dependence of the surface hydrogen abun- 
dance and the mass of the hydrogen envelope for the evo- 
lutionary stages following the end of the pre- WD evolu- 
tion is illustrated in Fig. 2] The role of element diffusion 

1 The occurrence of diffusion-induced hydrogen shell flashes is 
a theoretically known phenomena found in some detailed white 
dwarf evolutionary calculations (see for instance, Iben &; Mac- 
Donald 1986 in the context of intermediate-mass, carbon-oxygen 
white dwarfs). 
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Table 1. Selected stages for 0.33-Mq oxygen-core WD models considering element diffusion, 
as labeled in Fig.|3 From left to right: surface luminosity, effective temperature, age, surface 
gravity, nuclear luminosity, surface hydrogen abundance and hydrogen envelope mass. 





Log(L/L Q ) 


Log(T eff ) 


Age(10 6 yr) 


Log(g) 


Log(L nuc /L ) 


X H 


Log(M H /M ) 


(A) 


1.2715 


4.5640 


0.000 


5.8952 


2.843 


0.573 


-2.9003 


(B) 


-0.4906 


4.2337 


7.797007 


6.3360 


5.218 


1.000 


-3.4008 


(C) 


-1.4003 


3.8962 


7.797009 


5.8957 


4.657 


1.000 


-3.4008 


(D) 


3.4317 


3.7982 


7.797042 


0.6717 


3.2980 


0.304 


-3.4008 


(E) 


2.5554 


4.9168 


7.81 


6.0225 


2.544 


0.476 


-3.5674 


(F) 


0.0000 


4.5066 


7.98 


6.9369 


-1.499 


1.000 


-3.5987 


(G) 


-1.9967 


4.100 


184.2 


7.3095 


-2.968 


1.000 


-3.6502 


(H) 


-4.0289 


3.635 


4085 


7.4791 


-6.982 


1.000 


-3.6567 



Ages are counted from point A, which corresponds to an evolutionary stage just before the 
onset of the first thermonuclear flash. The surface hydrogen abundance is evaluated at a mass 
depth of 10~ 10 M„ below the stellar surface. 
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Figure 2. Hertzsprung-Russell diagram for the WD evolution of 
the 0. 33- Mq oxygen-core remnant when element diffusion is neg- 
elected. The WD experiences four hydrogen-shell flash episodes 
before reaching the final cooling branch. 



is clearly emphasized by this figure. During each short-lived 
flash episode, convection reaches deep enough into the star 
to dredge helium up, thus leading to envelopes made up by 
hydrogen and helium. This is responsible for the steep de- 
clines of the surface hydrogen abundance. For models with 
diffusion, after each flash, the purity of the outer layers is 
rapidly reestablished as soon as the star returns onto the 
cooling branch. Note that the time intervals during which 
the envelope remains helium-enriched are indeed extremely 
short. This is in contrast to the predictions of models with- 
out diffusion. For these diffusion-free models, there is no 
physical process other than convection that can transport 
hydrogen to the surface. Thus, the final hydrogen abundance 
(Xh ~ 0.35) is established by the mixing episode that oc- 
curred during the last flash. The mass of the hydrogen enve- 
lope is shown in the bottom panel of Fig. [I] An important 
feature revealed by this figure is that the remaining amount 
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Figure 3. Same as Fig.[2]but for the case when element diffusion 
is considered. Letters at various positions along the track indicate 
selected stages of the evolution described in Table 1. Also, the 
approximate location of the domain of ZZ Ceti instability strip is 
shown as a small shaded region. As a result of hydrogen chemically 
diffusing inwards during the cooling track, the WD experiences 
an additional thermonuclear flash, as compared with the situation 
of no diffusion. As a result of this diffusion-induced flash, the star 
for a brief time interval, is forced to evolve back to the red giant 
domain (point D). 



of hydrogen, shortly after hydrogen flashes have ceased, is 
markedly lower when diffusion is considered. Specifically, the 
hydrogen mass at the beginning of the final cooling track 
amounts to Mh = 2.52 x 1O _4 M0, as compared with the 
Mh = 3.82 x 10 -4 Mothat remains for the case when diffu- 
sion is neglected. The smaller hydrogen masses left in models 
with diffusion can be understood on the basis that, as the 
WD evolves along the cooling branch, the tail of the hy- 
drogen distribution chemically diffuses inwards where the 
temperature is high enough to burn it, a process which, as 
we mentioned, eventually induces the occurrence of an addi- 
tional thermonuclear flash. As can be seen, during this last 
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Diffusion 
No diffusion 



4 5 6 7 8 9 

Log(t WD ) [yr] 

Figure 4. Abundance by mass of hydrogen at a mass depth of 
10 — 10 M* below the stellar surface (upper panel) and mass of 
the hydrogen envelope in solar units (bottom panel) as a func- 
tion of WD age for the 0.33-MQOxygen WD model. Solid lines 
correspond to the case when element diffusion is considered and 
dashed lines when diffusion is neglected. At each flash episode, 
the surface abundance of hydrogen is abruptly reduced as a re- 
sult of convective mixing. Afterwards, diffusion causes the bulk of 
hydrogen to float on the surface again. Note also that during each 
flash episode, the mass of the hydrogen envelope is considerably 
reduced. Models with diffusion enter the final cooling branch (ages 
larger than 10 7 yr) with a markedly smaller hydrogen mass. By 
contrast models without diffusion retain a larger hydrogen en- 
velope which is largely burnt during the ensueing evolutionary 
stages. 



flash episode the total amount of hydrogen in the star is no- 
ticeably reduced as a result of the nuclear burning. As we 
shall see below, this fact produces a different cooling history 
even at very late stages of evolution. 

It is important to note that when diffusion is neglected, 
an appreciable amount of hydrogen is processed over the 
final cooling phase. Accordingly, we expect that the role 
played by stable hydrogen shell burning during the final WD 
cooling phase is different according to whether element dif- 
fusion is considered or not, an expectation that is indeed 
borne out by the results shown in Fig. |S] In the upper panel 
of this figure, we depict the ratio of nuclear-to-surface lumi- 
nosities as a function of WD age for the 0. 33- Mq oxygen-core 
WD. The solid line corresponds to the case when diffusion 
is included and the dashed line to the situation without dif- 
fusion. Only the evolutionary stages corresponding to the 
final cooling branch is depiced in the figure. Note that for 
these evolutionary stages of star models with diffusion, nu- 
clear energy release never constitutes the main contribution 
to the radiated luminosity. In particular, as the WD evolves 
across the domain of the ZZ Ceti stars, nuclear burning con- 
tributes at most 10 per cent to the surface luminosity. It is 
worth mentioning that the nuclear energy production for ef- 
fective temperatures lower than ~ 11 000K is almost entirely 
from the proton-proton chain. In fact, hydrogen burning via 
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l WD m 

Figure 5. Ratio of nuclear-to-surface luminosities (upper panel) 
and effective temperature (bottom panel) as a function of WD 
age for the 0.33-MQOxygen WD model. Solid lines correspond 
to the case when diffusion is considered and dashed lines when 
diffusion is neglected. Because the larger hydrogen envelope mass 
characterizing models without diffusion (see figure |3J, hydrogen 
burning for such models is the dominant source of energy even 
at advanced stages of evolution, with the result that evolution is 
considerably slowed down. 



CNO-cycle reactions abruptly ceases by the time the WD 
has approached the hot edge of the ZZ Ceti instability strip. 
Eventually, hydrogen burning becomes virtually extinct at 
the lowest luminosities that we computed. These results are 
in contrast to the situation in which element diffusion is 
neglected. Here, stable hydrogen burning via proton-proton 
chain becomes the dominant energy source even at advanced 
stages of evolution. Note that for this situation, hydrogen 
burning contributes almost 80-90 per cent to surface lumi- 
nosity even at T e g ~ 8000K. Thus, we conclude that element 
diffusion prevents hydrogen burning from being a main en- 
ergy source for most of the evolution of a low-mass, oxygen- 
core WD. This conclusion agrees with the predictions of Al- 
thaus et al. (2001ab) for the case of low-mass, helium-core 
WDs. 

Not surprisingly, the above-mentioned results have 
strong implications regarding the evolutionary ages of the 
models. This is documented in the bottom panel of Fig. [S] 
which illustrates the effective temperature versus the age of 
the WD. Again, the solid line corresponds to the evolution- 
ary sequence with element diffusion. Note that when diffu- 
sion processes are considered, cooling ages become roughly 
a factor 2-3 smaller. In fact, because nuclear reactions in 
models with diffusion do not dominate the energetics once 
the star has settled upon its final cooling track, to maintain 
its luminosity, the WD must extract energy from its relic 
thermal content, thus causing a much faster cooling. On the 
contrary, in the case of models without diffusion, evolution 
is dictated by residual hydrogen burning, giving rise to very 
long cooling ages. For instance, to reach T e g = 10 000 K, 
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Figure 6. Ratio of nuclear-to-surface luminosities (upper panel) 
and effective temperature (bottom panel) as a function of WD 
age for 0.33-MqWD sequences with diffusion. Solid (dashed) lines 
correspond to the case of oxygen-core (helium-core) models. For 
both cases, hydrogen burning do not represent the main source 
of energy. Note that oxygen-core WDs evolve much more rapidly 
than their helium-core counterparts. 



models with diffusion require 0.36 Gyr, whilst models with- 
out diffusion need about 1 Gyr. 

We compare now the evolutionary properties of our 
0.33-Afooxygen-core WD models with those of their helium- 
core counterparts. Specifically, the comparison is made with 
a 0.327Mq- helium WD sequence with diffusion as calcu- 
lated by Althaus et al. (2001a). We begin by examining the 
upper panel of Fig. HJ which illustrates the role of hydrogen 
burning as energy source for both sequences during the final 
cooling phase. For the stages depicted in the figure, nuclear 
burning is not the dominant source of star's luminosity for 
both core compositions (see Althaus et al. 2001a for the evo- 
lutionary properties of helium-core WDs when diffusion is 
considered); thus evolution is basically governed by the de- 
crease of the thermal heat content of the star. As is well 
known, for a given mass, an oxygen core has a lower total 
capacity of storing heat than a helium core (we remind the 
reader that the specific heat is proportional to the number of 
particles; hence, per gram of material, the specific heat is in- 
versely proportional to the atomic weight of the constitutive 
elements). Accordingly, we expect oxygen-core WDs to cool 
faster than helium-core ones. This expectation is confirmed 
by the bottom panel of Fig. HJ which shows the time taken 
by objects to cool down to a given efective temperature. We 
find that, at advanced stages, to reach a given temperature, 
the oxygen-core WD has to evolve in about half the time a 
helium-core WD need. 

The shape of the chemical composition profile is critical 
regarding the pulsational properties of WDs. In particular, 
it contributes to the shape of the Ledoux term appearing in 
the Brunt- Vaisala, frequency (Brassard et al. 1991) and plays 
a critical role in the phenomenon of mode trapping in WDs 



(see Brassard et al. 1992 and also Corsico et al. 2001). In Fig. 
|7|we show the abundance by mass of 1 H, 4 He, 12 C and ls O as 
a function of the outer mass fraction for the 0.33-A/ooxygen- 
core WD remnant. The upper panel depicts the abundance 
distribution shortly after the formation of the oxygen core 
and before the remnant experiences its first thermonuclear 
flash. The internal chemical profile emerging from the pre- 
WD evolution is characterized by an almost pure oxygen 
core of O.O8M0. The oxygen core is surrounded by a shell 
rich in helium, carbon, and oxygen with an overlying essen- 
tially pure helium buffer. The abundance distribution in the 
innermost regions is clearly different from that expected in 
the case of helium-core WDs, a feature which is expected to 
bear its signature in the theoretical period spectrum of these 
stars. The outermost layers are composed by helium and hy- 
drogen. The abundance of these elements is different from 
those assumed for the interstellar medium because mass loss 
during the close binary evolution exposed layers in which hy- 
drogen burning has occurred at earlier times. On the cooling 
track, the abundance distribution is altered by gravitational 
settling and chemical diffusion 2 . This is visualized in the 
bottom panel of Fig. [7| that illustrates the resulting chemi- 
cal profile when the WD has reached the domain of the ZZ 
Ceti instability strip on the final cooling branch. Note that 
element diffusion appreciably modifies the abundance distri- 
bution, particularly in the outer layers. Indeed, the hydrogen 
and helium profiles are markedly smoothed out by the time 
the ZZ Ceti domain is reached. At deeper layers, diffusion 
time-scale becomes much longer than the evolutionary time- 
scale and the chemical profile remains thus fixed during the 
entire WD evolution. 

The effect of diffusion on the element distribution 
within the star is emphasized in Fig. |H| where the abun- 
dances of 1 H, 4 He, 12 C, 14 N and ie O are shown as a func- 
tion of the outer mass fraction for O.33-M0Oxygen-core WD 
models at four selected epochs. The effect of diffusion on 
the element distribution is clearly noticeable. In particular, 
chemical diffusion moves a hydrogen tail inwards into hot- 
ter layers, causing the hydrogen burning rate to increase. 
At high T c ff values (panel b) this effect is responsible for 
the occurrence of the last flash episode. It is interesting to 
note that by the time the ZZ Ceti instability strip is reached 
(panel c) , there exists an appreciable tail of hydrogen in the 
helium buffer. In fact, by this time the hydrogen distribution 
reaches its maximum depth. At the same time, there are no 
metals in the outermost 2 x 10~ 4 M0of the star due to the 
occurrence of gravitational settling. 



4 PULSATIONAL ANALYSIS 

We have seen that the evolutionary properties of the oxygen- 
core WDs differ appreciably from those of helium-core WDs, 
a fact which could eventually be used to shed light on the 
internal composition of observed low-mass WDs. However, 
a more promising and potential way of distinguishing both 



2 It should be kept in mind that an appreciable fraction of the 
hydrogen envelope is consumed by nuclear burning during the 
flash episodes. 
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Figure 7. Chemical abundance distribution in terms of the outer 
mass fraction for the 0.33-AfQOxygen-core WD. Upper panel cor- 
responds to an evolutionary stage shortly after the formation of 
the oxygen core and bottom panel to an evolutionary stage near 
the ZZ Ceti instability strip. Note that element diffusion markedly 
modifies the shape of the hydrogen and helium profiles. 



types of stars would be by means of the study of their pulsa- 
tional patterns if low-mass WDs were in fact found to pul- 
sate as ZZ-Ceti type variables. Although signs of variability 
in low-mass WDs have not been detected so far, we judge 
that an exploration of their oscillatory properties is never- 
theless worthwhile. Indeed, in the context of more massive 
WDs, it is well known that asteroseismology is a very pow- 
erful tool for probing into the internal composition of de- 
generate stars (see, e.g., Bradley 1998 and Metcalfe 2003). 
In particular, we are interested here in studying both the 
adiabatic and non-adiabatic pulsational properties. 



4.1 Adiabatic pulsations 

In order to explore the adiabatic pulsational properties of 
the WD models presented in this work, we have employed 
a general Newton-Raphson pulsational code coupled to the 
LPCODE evolutionary code as detailed in Corsico et al. 
(2001, 2002; see also Corsico 2003). In particular, the com- 
putation of the Brunt- Vaisala frequency (N), an issue of 
fundamental relevance in the context of pulsating WD, is 
performed by adopting the prescription suggested by Bras- 
sard et al. (1991). Specifically, our code computes the so 
called "Ledoux term" (B), given by 



B = - 


1 

Xt 


with 






/ dlnP 




VainAi 


and 





din A; 



(3) 



(4) 



Any local feature in the profiles of the internal chemi- 
cal abundances is reflected by the Brunt- Vaisala frequency 
through the Ledoux term, as it can be clearly noted from 
the following expression for N (Brassard et al. 1991): 



N 2^9_^Xr (Vad — V + -B) . 
" Xp 



(6) 



Here, Xi is the mass fraction of specie i, n the number of 
ionic species present in the plasma, g the local gravity, V a d 
and V the adiabatic and actual temperature gradients, re- 
spectively, and 



dlnP 
d\np 



(7) 



T,{Xi} 



In this section we explore the adiabatic pulsational 
properties of low-mass, oxygen-core WD models and we 
compare them with those corresponding to the helium-core 
counterparts. Specifically, we have computed the adabatic 
pulsational pattern of gravity modes corresponding to an 
evolutionary sequence of 0.33- M@ oxygen-core WD models 
for an effective temperature range of 15 000-8000 K. Addi- 
tionally, we have carried out pulsational calculations cor- 
responding to an evolutionary sequence of 0. 33- Mq helium- 
core WD models for the same T e g interval. To assess the 
mode-trapping/confinement features associated with both 
kind of objects, we restricted ourselves to analyzing two 
selected models picked out from both sequences at T e g « 
10 000 K. Hereafter, we shall denote these template mod- 
els as OCWD (oxygen-core WD) and HECWD (helium- 
core WD) models. We begin by examining the shape of the 
Ledoux term and the Brunt- Vaisala frequency in Figs.[§]and 
1101 f middle and bottom panels, respectively). To make eas- 
ier the association amongst the various chemical transitions 
and the corresponding features in N, we have also included 
in these figures the internal chemical profiles (upper panels). 

For the case of the OCWD model, the most outstand- 
ing characteristic in the shape of N is the peaked feature 
at log(l - M r /M,) » -0.2 (Fig. EJ. This is directly asso- 
ciated with the abrupt variation in the oxygen abundance 
at such very deep layers. On the other hand, the oxygen- 
helium interface is responsible for the less narrow neighbour- 
ing peak in N, and finally, at log(l — M r /M,) ~ —3.25 there 
is a broad, smooth bump which is caused by the helium- 
hydrogen transition region. The shape of N correspond- 
ing to the HECWD model is documented in the bottom 
panel of Fig. 1101 As can be appreciated from the figure, in 
this case the only local feature present in N is that result- 
ing from the helium-hydrogen chemical transition zone. The 
comparison of the Brunt- Vaisala frequency for both mod- 
els deserves some comments. As is well known, any local 
feature in the Brunt- Vaisala frequency constitutes a source 
for trapping or confinement of modes (see Brassard et al. 
1992; Corsico et al 2002). As a result, we expect the mode- 
trapping/confinement properties of OCWD model to be con- 
siderably more featured than those of the HECWD model. 
The other point worthy of comment is the fact that N in the 
HECWD model shows slightly greater values at the inner- 
most core region than for the OCWD case. Although this 
seems to be a subtle difference, it causes the periods to be 
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Figure 8. Chemical abundance distribution for 1 H, 4 He, 12 C, 14 N and le O as a function of the outer mass fraction for the 0.33- 
M0Oxygen-core WD at four selected evolutionary stages. Panel a shows the chemical stratification during the second flash episode, 
shortly after the occurrence of the outer convective mixing. Panel b corresponds to the stage just before the occurrence of the last flash. 
To emphasize the importance of chemical diffusion, we also show with thin solid line the hydrogen profile emerging from the previous 
flash. Panels c and d correspond to models on the final cooling branch: the ZZ Ccti domain and the last computed model, respectively. 
Note the role of both gravitational settling and chemical diffusion. In particular, chemical diffusion causes some hydrogen to reach deeper 
layers. At high T c g values (panel b) this effect is responsible for the occurrence of the last flash episode. Values of (log(L/Lo) and 
logT eff ) for panels a, b, c and d are, respectively: (0.96; 4.34), (0.72; 4.64), (-2.53; 3.98) and (-4.03; 3.63). 



shorter and the asymptotic period spacing to be somewhat 
smaller in the case of the OCWD model. 

In what follows we comment on the results we obtained 
for the pulsational properties of our models. Specifically we 
have computed the adiabatic, non-radial g-modes of low 
degree (£ = 1,2,3) with overtones covering a wide period 
window of 50 — 2000 sec. In addition, related quantities 
such as oscillation kinetic energy (Bkin) and first-order rota- 
tional splitting coefficient (CV,fc) have been obtained for each 
computed mode. We mention that the radial eigenfunction 
(Sr/r) is normalized to 1 at the surface. The main observable 
quantities in pulsating WD stars are the periods (Pk) and 
the forward period spacing (APk = Pk+i — Pk)- Thus, we 



shall employ these quantities as the main discriminants to 
demonstrate the different pulsational behaviours associated 
with OCWD and HECWD model stars. 

We begin by examining Fig. II II which shows the periods 
as a function of the overtone number for £ —1,2 and 3 corre- 
sponding to OCWD (thick lines) and HECWD (thin lines) 
models. Note that all the periods corresponding to HECWD 
model are longer than for the case of the OCWD one. As 
mentioned, this is a consequence of N being slightly larger 
at the innermost core region in the case of the HECWD 
model. Note also for this model the near uniformity in the 
period distribution as compared with the OCWD case (see 
below) . 
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Figure 9. The internal chemical profile, the Ledoux term and 
the Brunt- Vaisala frequency (upper, middle and bottom panel, 
respectively) in terms of the outer mass fraction, corresponding 
to a selected O.33-M oxygen-core WD model at T cff rs 10000 K. 
Note that N exhibits a peaked feature at the innermost region of 
the model. It is a direct consequence of the steep change of the 
16 abundance at log(l — Mr/M,) ss -0.2. The remainder, less 
narrow features, are associated with the oxygen-helium and the 
helium-hydrogen interfaces. 
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Figure 10. Same as Fig. [5] but for a helium-core WD model of 
0.33-AfQat the same effective temperature. In this case, the only 
feature characterizing the shape of the Brunt- Vaisala frequency 
is that associated with the helium-hydrogen interface (the only 
chemical transition region in this model). 
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Figure 11. Periods as a function of the overtone number for I 
=1, 2 and 3 corresponding to OCWD (thick lines, fill dots) and 
HECWD (thin lines, hollow dots) models. 



In Fig. El w e show APk, log(_Ekin) and Ci t k (upper, 
middle and bottom panels, respectively) in terms of peri- 
ods, for 1= 1,2,3 (left, centre and right) corresponding to 
the OCWD model. Fig. 1131 depicts the same quantities but 
for the case of the HECWD model. We begin by examining 
the APk distribution in Fig. 1121 As expected, the spacing 
of consecutive periods is markedly featured, reflecting the 
three sources of trapping/confinement of modes correspond- 
ing to each of the three chemical interface. An inspection of 
the figures reveals that the amplitude of APk corresponding 
to the OCWD model is typically larger as compared with 
the HECWD model. Note also the strong minima in APk 
in the OCWD model. These minima correspond to modes 
partially confined to the deepest regions of the star. Indeed, 
at the high-density zone bounded by the stellar centre and 
log(l — M r /M*) ~ —0.2, such modes have larger amplitudes 
of their eigenfunctions as compared with those of the neigh- 
boring overtones. Consequently, such modes must exhibit 
local maxima in the -Ekin distribution (see middle panels of 
Fig. I12H . In contrast, the much less pronounced minima ex- 
hibited by APk in the HECWD model are related to trapped 
modes in the outer layers. Thus, these modes are character- 
ized by low kinetic energy values (Fig. 1131 . Finally, we note 
that the distribution of the Ce t k values corresponding to the 
HECWD model exhibits less structure as compared with the 
case of the OCWD one, particularly regarding to the short- 
and intermediate-period domains. 

Clearly, from the above considerations we conclude that 
there exist marked differences beetwen the adiabatic pulsa- 
tional properties of helium- and oxygen-core low-mass WDs. 
This in principle turns the asteroseismological techniques 
usually employed in variable WD studies into an important 
tool for distinguishing both types of stars. 
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Figure 12. The forward period spacing (upper panels), the logarithm of the kinetic energy (middle panels) and the first order rotational 
splitting coeficients (lower panels), in terms of the periods for (/-modes with 1 = 1 (left), 1 = 2 (centre) and t = 3 (right), for the same 
WD model analyzed in Fig. [5] (OCWD models). In the interest of clarity we have ommited the symbols representing pulsation modes. 
Dotted lines represent the asymptotic period spacing (upper panels) and the asymptotic value of Ci k (lower panels). 



4.2 Non-adiabatic pulsations 

Finally, we comment on the main results from non-adiabatic 
pulsation calculations of the oxygen-core WD models. We 
emphasize that for the WD regime, we modelled convec- 
tion according to the mixing-length theory setting the mix- 
ing length to 1.5 pressure scale-heights. The grey domain 
in Fig. 1141 shows the extension of the superficial convection 
zone as a function of effective temperature. Note that at 
T e ff « 10 000 K, below which modes are found unstable (see 
below), the base of the convection zone sinks appreciably 
into the star. The linear non-adiabatic oscillation spectra 
were computed for dipolar (I = 1) p-modes in the period 
range between 100 and about 1200 sec using the same Ric- 
cati shooting- method as referred to in Gautschy et al. (1996) 
(see also Gautschy & Althaus 2002). The non-adiabatic com- 
putations were done assuming a frozen-in convective flux. 



Figure 1151 displays the modal diagram for the 0.33- 
M©oxygen-core WD sequence. Specifically, the figure de- 
picts the periods of the lowest order £ = 1 modes as a 
function of effective temperature. Overstable modes are in- 
dicated by black dots. Like radial orders are connected by a 
continuous line. The highest order shown is k — 18. An in- 
stability domain is found that starts at log Teg = 3.993 and 
continues below the cool end of the model series analyzed. 
Only modes with periods exceeding about 680 s become un- 
stable. A striking feature concerns the position of the blue 
edge: It does not depend on the radial order, i.e. on the 
length of the periods. A survey computation showed further- 
more that the saturation of the excitation must be beyond 
k = 30. In other words, we did not find the long-period edge 
of the instability domain. 

Figure Hoi shows the imaginary parts of the eigenvalues 
as a function of T e a. Note that the maximum of instabil- 
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Figure 13. Same as Fig. I12l but for the model analyzed in Fig llUI (HECWD model). 



ity, i.e. the minimum of o\, is encountered very close to the 
blue edge. The excitation rates decline quickly toward lower 
temperatures. Notice, as referred to in the last paragraph, 
that the instability (i.e., the magnitude of <j\) continues to 
get stronger to higher overtone numbers. 

We mention that the non-adiabatic computations in- 
cluded the effect of nuclear burning, i.e. the e-mechanism. 
However, we considered only the equilibrium situation, ne- 
glecting the perturbation of the abundance of the reactants. 
Wc find that the work integrals have no appreciable contri- 
bution from e-mechanism. A point of interest is the question 
if pulsations might trigger a very late hydrogen shell flash. 
Based on the fact that the eigenfunctions have very low am- 
plitudes in hydrogen-burning regions, we do not expect a 
pulsation-induced thermonuclear flash to be feasible 3 . 



3 The g-mode pulsations' preponderant horizontal displacements 
are unlikely to disturb the geometry of the nondegenerate H- 
burning shell such as to trigger a thermal runaway. 



Finally, since no convection - pulsation interaction was 
considered, the non-adiabatic results indicate at the very 
best that the k mechanism destabilizes low-degree modes of 
intermediate order. The action of convection can still modify 
the picture considerably, even on a qualitative level. 



5 CONCLUSIONS 

Motivated by theoretical evidence suggesting that some of 
the presumed low-mass, helium-core WDs could actually be 
WDs with oxygen cores (Iben & Tutukov 1985 and Han, 
Tout & Eggleton 2000) , we explored in this paper the evolu- 
tion of low-mass WDs with oxygen cores in a self-consistent 
way with the predictions of nuclear burning, the chemical 
evolution caused by time-dependent element diffusion, and 
the history of the white dwarf progenitor. We concentrated 
on an oxygen-core WD remnant of mass 0.33 M@. The evo- 
lutionary stages prior to the white dwarf formation were 
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Figure 14. Location of the convection zone as given by the outer 
mass fraction as a function of effective temperature (grey do- 
main). The vertical, dashed line indicates the value of effective 
temperature below which g-modes become unstable. The thick 
line marks the outer edge of the convection zone which is near 
the location of the photosphere. 
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Figure 15. Modal diagram for I = 1 of the 0.33-MQOxygen- 
core white dwarf sequence. The dots on some intermediate-order 
modes indicate their instability. 
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Figure 16. Same of Figure 1151 but for the imaginary part of 
selected eigenfrequencies. Negative numbers mean oscillatory in- 
stability. The radial order of the eigenmodes is labelled on the 
graph. 

also fully accounted for by computing the conservative bi- 
nary evolution of an initially 2.5-MePop. I star with a 1.25 
MQCompanion, and an orbital period of Pi= 3 days. In ad- 
dition, we explored the adiabatic and non-adiabatic pulsa- 
tional properties of the resulting WD models. 

As for the main evolutionary results we mention: 

- The star's evolution is characterized by the occurrence 
of several episodes of thermonuclear flashes due to unstable 
hydrogen burning at the base of the hydrogen-rich envelope. 

- Diffusion processes play an important role in the evo- 
lution of low-mass oxygen-core WDs. In particular, the re- 
maining amount of hydrogen, shortly after hydrogen flashes 
have ceased, is markedly lower when diffusion is considered. 
This fact produces a different cooling history even at very 
late stages of the evolution. In particular, element diffu- 
sion prevents hydrogen burning from being the main energy 
source for most of the evolution of a low-mass, oxygen-core 
WD. As a result, cooling ages are roughly a factor 2 — 3 
smaller than ages derived from models neglecting diffusion. 

- As compared with the helium-core counterparts, low- 
mass, oxygen-core WDs cool much faster once they have 
settled on their final cooling track. In particular, at advanced 
stages, to reach a given temperature, the oxygen-core WD 
evolves about twice as fast as a helium-core WD. 

The main pulsational properties of our models are: 

- The structure of the pulsational spectrum of the 
oxygen-core WD models is notably richer than that of the 
helium-core ones. This fact is apparent from the period 
spacing diagrams, which exhibit large amplitudes of trap- 
ping/confinement imprints. 

- The presence of an oxygen core is responsible for some 
partially confined modes that are characterized by compar- 
atively large values of kinetic energy and strong minima in 
a period-spacing diagram. 

- Modes with periods longer than about 680 sec become 
unstable at effective temperatures below about 10000 K. 

We conclude that both the pulsational and evolution- 
ary properties of low-mass, carbon-oxygen WDs turn out to 
be sustantially different from those of the helium-core WDs. 
In particular, our results suggest that a discrimination be- 
tween helium-core and carbon/oxygen-core WDs by means 
of asteroseismology is possible if low-mass WDs were in fact 
found to pulsate as ZZ Ceti - type variables. Despite the lack 
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of signs of variability in low-mass WDs so far, we consider 
the monitoring of low-mass WDs to find variable candidates 
to be worthwhile. 

Complete tabulations of the results of the present paper 
are available at http://www.fcaglp.unlp.edu.ar/evolgroup/ 
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support by the Instituto de Astroffsica de La Plata. This re- 
search has made use of NASA's Astrophysical Data System 
Abstract Service. 
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